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^ ■ Abstract 
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T— I , We examine the use of two superbeam neutrino oscillation experiments 

\ with baselines < 1000 km to resolve parameter degeneracies inherent in the 

' three-neutrino analysis of such experiments. We find that with appropriate 

Q_i' choices of neutrino energies and baselines two experiments with different base- 

^ • lines can provide a much better determination of the neutrino mass ordering 

D i than a single experiment alone. Two baselines are especially beneficial when 

\ the mass scale for solar neutrino oscillations ^m^^^ is > 5 x 10~^ eV^. We also 

^ ^ ' examine CP violation sensitivity and the resolution of other parameter degen- 

^ ■ eracies. We find that the combined data of superbeam experiments with base- 

^ . lines of 295 and 900 km can provide sensitivity to both the neutrino mass or- 

dering and CP violation for sin 2013 down to 0.03 for \5ml^J\ ~ 3x lO'^ eV^. 
It would be advantageous to have a 10% determination of |<5?n^^;jj| before the 
beam energies and baselines are finalized, although if is not that well 

known, the neutrino energies and baselines can be chosen to give fairly good 
sensitivity for a range of |(5m^^jn|. 



I. INTRODUCTION 



Atmospheric neutrino data from Super-Kamiokande provides strong evidence that i/^'s 
created in the atmosphere oscillate to 1/^ with mass-squared diflFerence ~ 3 x 10~^ eV^ 

and almost maximal amplitude [1]. Furthermore, the recent solar neutrino data from the 
Sudbury Neutrino Observatory (SNO) establishes that electron neutrinos change flavor as 
they travel from the Sun to the Earth: the neutral-current measurement is consistent with 
the solar neutrino flux predicted in the Standard Solar Model [2], while the charged-current 
measurement shows a depletion of the electron neutrino component relative to the total 
flux [3]. Global fits to solar neutrino data give a strong preference for the Large Mixing 
Angle (LMA) solution to the solar neutrino puzzle, with 5m1^] ~ 5 x 10~^ eV^ and amplitude 
close to 0.8 [3,4]. 

The combined atmospheric and solar data may be explained by oscillations of three 
neutrinos, that are described by two mass-squared differences, three mixing angles and a 
CP violating phase. The atmospheric and solar data roughly determine 5rnl^^, <^"^soi ^^'^ 
corresponding mixing angles. The LMA solar solution will be tested decisively (and Sml^Y 
measured accurately) by the KamLAND reactor neutrino experiment [5,6]. More precise 
measurements of the other oscillation parameters may be performed in long-baseline neutrino 
experiments. The low energy beam at MINOS [7] plus experiments with ICARUS [8] and 
OPERA [9] will allow an accurate determination of the atmospheric neutrino parameters 
and may provide the first evidence for oscillations of — > v^. at the atmospheric mass 
scale [10]. It will take a new generation of long-baseline experiments to further probe — > Ve 
appearance and to measure the leptonic CP phase. Matter effects are the only means to 
determine sgn{5m\^^] once sgn((5m^j^) is known, the level of intrinsic CP violation may 
be measured. Matter effects and intrinsic CP violation both vanish in the hmit that the 
mixing angle responsible for — > oscillations of atmospheric neutrinos is zero. 

It is now well-known that there are three two-fold parameter degeneracies that can occur 
in the measurement of the oscillation amplitude for z/^ appearance, the ordering of the 

neutrino masses, and the CP phase [11]. With only one v and one v measurement, these 
degeneracies can lead to eight possible solutions for the oscillation parameters; in most cases, 
CP violating {CPV) and CP conserving (CPC) solutions can equally explain the same data. 
Studies have been done on how a superbeam [11-16], neutrino factory [16-18], superbeam 
plus neutrino factory [19], or two supcrbcams with one at a very long baseline [20,21] could 
be used to resolve one or more of these ambiguities. 

In this paper we show that by combining the results of two superbeam experiments 
with different medium baselines, < 1000 km, the ambiguity associated with the sign of 
^^'^atm '^^^ t)e resolved, even when it cannot be resolved by the two experiments taken 
separately. Furthermore, the ability to determine sgn(5m^^.jjj) from the combined data is 
found to not be greatly sensitive to the size of Sm^^^i, unlike the situation where data from 
only a single baseline is used. If both experiments are at or near the peak of the oscillation, 
a good compromise is obtained between the sensitivities for resolving sgn(5m^^.j^) and for 
establishing the existence of CP violation. If l^^atml ^'^^ known accurately, the neutrino 
energies and baselines can be chosen to give fairly good sensitivity to the sign of Sm^^^^ and 
to CP violation for a range of |5/«atml- 

The organization of our paper is as follows. In Sec. II we discuss the parameter de- 
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generacies that can occur in the analysis of long-basehne oscillation data. In Sec. Ill we 
analyze how two long-baselinc siipcrbeam experiments can break degeneracies, determine 
the neutrino parameters, and establish the existence of CP violation in the neutrino sector, 
if it is present. A summary is presented in Sec. IV. 

II. PARAMETER DEGENERACIES 

We work in the three-neutrino scenario using the parametrization for the neutrino mixing 
matrix of Ref. [11]. If we assume that z/3 is the neutrino eigenstate that is separated from 
the other two, then = ^w^atm ^"^^ ^^e sign of Sm}^^ can be either positive or negative, 

corresponding to the mass of 2/3 being either larger or smaller, respectively, than the other 
two masses. The solar oscillations are regulated by 5m|i = ^f^l^i, and thus |5m|i| <S 
If we accept the likely conclusion that the solar solution is LMA [3,4], then Sm^i > and we 
can restrict 612 to the range [0, vr/4]. It is known from reactor neutrino data that ^13 is small, 
with sin^ 2^13 < 0.1 at the 95% C.L. [22]. Thus a set of parameters that unambiguously 
spans the space is Sml^ (magnitude and sign), Sm^i, sin^ 26*12, sin 6*23, and sin^ 26*13; only the 
023 angle can be below or above 7r/4. 

For the oscillation probabilities for — > Ue a-nd z/^ Pe we use approximate expressions 
given in Ref. [11], in which the probabilities are expanded in terms of the small parameters ^13 
and 617121 [23,24], which reproduces well the exact oscillation probabilities for > 0.5 GeV, 
^13 ^ 9°; and L < 4000 km [11]. In all of our calculations we use the average electron 
density along the neutrino path, assuming the Prehminary Reference Earth Model [25]. Our 
calculational methods are described in Ref. [12]. 

We expect that |5?Ti|i| and sin^ 2^^23 will be measured to an accuracy of 10% at 3a 
from survival in long-baseline experiments [7-10], while 5m\i will be measured to 

an accuracy of ~ 10% at 2a and sin^ 26*12 will be measured to an accuracy of ±0.1 at 2a in 
experiments with reactor neutrinos [6]. The remaining parameters (^13, the CP phase 5, and 
the sign of 5m|i) must be determined from long-baseline appearance experiments, principally 
using the modes z/g and z/^ z/g with conventional neutrino beams, or z/g z/^ and 

i^fj. at neutrino factories. However, there are three parameter degeneracies that can 
occur in such an analysis: (i) the {6,613) ambiguity [17], (ii) the sgn(5m|i) ambiguity [13], 
and (iii) the (6*23, 7r/2 — ^23) ambiguity [11] (see Ref. [11] for a complete discussion of these 
three parameter degeneracies) . In each degeneracy, two different sets of values for 5 and ^13 
can give the same measured rates for both u and P appearance and disappearance. For each 
type of degeneracy the values of ^^13 for the two equivalent solutions can be quite different, 
and the two values of 6 may have different CP properties, e.g., one can be CP conserving 
and the other CP violating. 

A judicious choice of L and E^, can reduce the impact of the degeneracies. For example, 
if L/E^ is chosen such that A = [5m|i|L/(4£^j,) = 7r/2 (the peak of the oscillation in 
vacuum), then the cos 5 terms in the average appearance probabilities vanish, even after 
matter effects are included [11]. Then since it is sin 5 that is being measured, the (5, 6*13) 
ambiguity is reduced to a simple {S,tt — S) ambiguity, CPV solutions are no longer mixed 
with CPC solutions, and ^13 is in principle determined (for a given sgn(5m3^) and ^23)- 
If L is chosen to be long enough (> 1000 km), then the predictions for 5m|i > and 
Smli < no longer overlap if ^13 > a few degrees, and the sgn((5m|i) ambiguity is removed; 
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our previous studies indicated that for ^mli = 5 x 10^'' eV^ this happens at L > 1300 km 
if sin^ 2^^i3 > 0.01 [11] (before experimental uncertainties are considered). However, the 
persistence of the sgn((5m|i) ambiguity is highly dependent on the size of the solar oscillation 
mass scale, because large values of Sml^ cause the predictions for Sml^ > and Sml^ < 
to overlap much more severely than when 5mli is smaller. Also, existing neutrino baselines 
are no longer than 735 km. In this paper we explore the possibility that two experiments 
with medium baselines (< 1000 km) can determine sgn(5m|^), even when data from one of 
the baselines alone cannot. We then address the sensitivity for establishing CP violation. 



III. JOINT ANALYSIS OF TWO SUPERBEAM EXPERIMENTS 
A. Description of the experiments and method 

For our analysis we take one baseline to be 295 km, the distance for the proposed experi- 
ment from the Japan Hadron facility (JHF) to the Super-Kamiokandc detector at Kamioka. 
For the neutrino spectrum of this experiment we use their 2° off-axis beam with average 
neutrino energy of 0.7 GeV [26]. For the second experiment we assume an off-axis beam in 
which the beam axis points at a site 735 km from the source (appropriate for a beamline 
from NuMl at Fermilab to Soudan, or from CERN to Gran Sasso). For the off- axis spectra 
of the NuMl experiment we use the results presented in Ref. [27], which provides neutrino 
spectra for 39 different off-axis angles ranging from 0.32° to 1.76°. 

Using the off-axis components of the beam has the advantage of a lower back- 
ground [15,28,29] due to reduced Ug contamination and a smaller high-energy tail. Off-axis 
beams also offer flexibility in the choice of L and E,^. For example, for a beam nominally 
aimed at a ground-level site a distance Lq from the source, the distance to a ground-level 
detector with off-axis angle 6oa can lie anywhere in the range 

2Re sm{e - doA) <L<2Re sm{e + doA) , (1) 

where sin e = Lo/(2i?e), and Re = 6371 km is the radius of the Earth. Then for < Rj 
the possible range of distances for an off-axis detector at approximately ground level is 

Lq — 2Re9oA ^ L < Lo + 2Re9oA ■ (2) 
The neutrino energy and neutrino flux decrease with increasing off-axis angle as 

0.43K ^ E?; 

where 7 = E^/m.,^ is boost factor of the decaying pion. Thus a wide range of L and E^, 
can be achieved with a single flxed beam, although the event rate will drop with increasing 
off-axis angle because the flux decreases and the neutrino cross section is smaller at smaller 
El, (thereby putting a limit on the usable range of L and E^) . 

For the first experiment at Li = 295 km, we assume that the neutrino spectrum is chosen 
so that the cos 6 terms in the u and u oscillation probabilities vanish (after averaging over the 
neutrino spectrum), using the best existing experimental value for 5ml^. The JHF 2° off-axis 
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beam [30] satisfies this condition for Sm^i = 3 x 10"'^ eV^. This spectrum choice reduces the 
{S, 9iQ,) ambiguity to a simple (5, tt — 5) ambiguity, as described in Sec. II. For the second 
experiment we allow L2 and Oqa to vary within the restrictions of Eq. (2). This flexibility can 
be fully utilized if a deep underground site is not required; the short duration of the beam 
operation (an 8.6 /is pulse with a 1.9 s cycle time [31]) may enable a sufficient reduction 
in the cosmic ray neutrino background. We assume that the proton drivers at the neutrino 
sources have been upgraded from their initial designs (from 0.8 to 4.0 MW for JHF [30] and 
from 0.4 to 1.6 MW for FNAL [32]), so that they are both true neutrino superbeams. We 
assume two years running with neutrinos and six years with antineutrinos at JHF, and two 
years with neutrinos and five years with antineutrinos at FNAL; these running times give 
approximately equivalent numbers of charged-current events for neutrinos and antineutrinos 
at the two facilities, in the absence of oscillations. For detectors, we assume a 22.5 kt 
detector in the JHF beam (such as the current Super-K detector) and a 20 kt detector 
in the FNAL beam (which was proposed in Ref. [15]). Larger detectors such as Hyper- 
Kamiokande or UNO would allow shorter beam exposures or higher precision studies. In all 
of our calculations, we assume \Sml^\ = 3 x 10"^ eV^, 623 — 7r/4, Sml^ — 5 x 10"^ eV^, and 
sin^ 29i2 = 0.8, unless noted otherwise. 

We first consider the minimum value of sin^ 2^13 for which the signal in the neutrino 
appearance channel can be seen above background at the 3a level (the discovery reach), 
varying over a range of allowed values for Oqa and L2 in the second experiment. The 
discovery reach depends on the value of S and the sign of 5m|^; the best (when Sm^i > 0) 
and worst (when Sm'^i < 0) cases in the u channel (after varying over S) are shown in 
Fig. 1. In the u channel, the best case occurs for Sni^i < and the worst for > 0. In 

our calculations we assume a background that is 0.5% of the unoscillated charged-current 
rate (see Ref. [15]), and that the systematic error is 5% of the background. However, we 
note that our general conclusions are not significantly affected by reasonable changes in 
these experimental uncertainty assumptions. Detector positions where there is no cos 5 
dependence in the rates arc denoted by boxes. The best reach is sin^ 29ri ~ 0.003, which 
occurs for Oqa — 0.5-0.9°. In the worst case scenario the reach degrades to sin^ 29i2 — 0.01. 

The measurement of P and P at Li allows a determination of sin^ 2^13 and sin 5, modulo 
the possible uncertainty caused by the sign of f^n?!^, assuming for the moment that ^23 = 7''/4, 
so there is no (^^23,7r/2 — ^^23) ambiguity. The question we next consider is whether an 
additional measurement of P and P at L2 can determine sgn(5m|^), measure CP violation, 
and distinguish 6 from n — 6. We define the of neutrino parameters {6', 9[^) relative to 
the parameters {S, ^13) as 

where Ni and N- are the event rates for the parameters (5, ^13) and {5',9[^), respectively, 
6Ni is the uncertainty in Ni, and i is summed over the measurements being used in the 
analysis {u and u at Li and u and u at L2). For SNi we assume that the statistical error 
for the signal plus background can be added in quadrature with the systematic error. For a 
two-parameter system {S and ^13 unknown), two sets of parameters can be resolved at the 
2(7 (3(7) level if > 6.17 (11.83). 
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B. Determining the sign of Srn^i 

To determine if measurements at Li and L2 can distinguish one set of oscillation param- 
eters with one sign of Sml^ from all other possible sets of oscillation parameters with the 
opposite sign of Sm^i, we sample the (5, ^13) space for the opposite sgn(5m|i) using a fine 
grid with 1° spacing in 5 and approximately 2% increments in sin^ 2^13. If the between the 
original set of oscillation parameters and all of those with the opposite sgn(5m3^) is greater 
than 6.17 (11.83), then sgn{6ml^) is distinguished at the 2a (Scr) level for that parameter 
set. 

Figure 2 shows contours (in the space of possible L2 and 9oa for the second experiment) 
for the minimum value of sin^ 2^13 (the sin^ 2^13 reach) for distinguishing sgn{Sm^{) at the 
3(T level when u and u data from Li and L2 are combined. As in Fig. 1, the boxes indicate 
the detector positions where the cos 5 terms in the average probabilities vanish. The best 
reach of about sin^2^^i3 ~ 0.03 can be reahzed for Oqa — 0.7-1.0° and L2 values near the 
maximum allowed by Eq. 2 (~ 875-950 km). Table I shows the sensitivity for determining 
sgn(5m|]^) for different combinations of detector size and proton driver power in the two 
experiments. The table shows that once enough statistics are obtained at JHF (with a 
22.5 kt detector and a 4 MW source), combined JHF and NuMI data significantly improve 
the sin^2^i3 reach for determining sgn(5m3j^) at 3cr (by nearly a factor of two compared to 
data from a 1.6 MW NuMI alone). 



TABLE I. sin^ 2^13 reach for determining the sign of Sm^i at 3a using u and u data from JHF 
at 295 km and NuMI at L2, for various detector sizes and proton driver powers. The approximate 
range of Oqa that can obtain the reach shown is given in parentheses; L2 ~ 900 km in all cases. 







NuMI (20 kt) 




JHF 


0.4 MW 




1.6 MW 




sin2 2013 {Ooa) 




sin2 2^?i3 (eoA) 


no JHF data 


0.09 (0.7-1.0°) 




0.05 (0.8-1.0°) 


22.5 kt, 0.8 MW 


0.07 (0.8-1.0°) 




0.04 (0.9-1.0°) 


22.5 kt, 4.0 MW 


0.06 (0.7-1.0°) 




0.03 (0.7-1.0°) 


450 kt, 4.0 MW 


0.05 (0.6-1.0°) 




0.02 (0.7-0.9°) 



The ability to distinguish the sign of Smli is greatly affected by the size of the solar 
mass scale Smli, because the predictions for Smli > and Smli < overlap more for larger 
values of Smli- In Fig. 3a we show the region in (5, sin^2^i3) space for which parameters 
with 5mli > can be distinguished from all parameters with 5ml^ < at the 3a level 
for several possible values of Srri^i, using combined data from Li = 295 km and L2 = 
890 km, with Oqa = 0.74° for the second experiment. With this configuration the cos 5 
terms in the average probabilities vanish for both experiments and nearly maximal reach for 
distinguishing sgn{6mli) is achieved. A similar plot using only data at L2 = 890 km and 
9oA — 0.74° is shown in Fig. 3b. We do not show a corresponding plot for Li — 295 km 
because the shorter baseline severely inhibits the determination of sgn((5m|i). A comparison 
of the two figures shows that for 5 — 270° (where the Sml^ > predictions have the 



6 



least overlap with any of those for Sm^i < 0) the sensitivity to sgn(5m|]^) is not significantly 
improved by adding the data at Li. However, at 6 = 90° the ability to distinguish sgn(5m|;^) 
is much less affected by the value of Sml^ when the data at Li is included. With data only 
at L2, sgn(5m3^) can be determined for sin^ 2^^i3 = 0.1 when 6 — 90° only for Sm^^ < 
8 X 10~^eV^, while with data at Li and L2 it can be determined for sin^ 2^13 as low as 0.04 
for Sm2i as high as 2 x 10~^ eV^. The corresponding results for 5m^i < are approximately 
given by reflecting the curves in Fig. 3 about d = 180°. 

We conclude that combining measurements of — > i^g and — > from two superbeam 
experiments at different L results in a much more sensitive test of the sign of Smli than one 
experiment alone, especially for larger values of the solar mass scale (5m|^ . 

The ability to determine sgn((5m|]^) is also affected by the value of 623. We found that the 
sin^ 2^13 reach for determining sgn(5m|]^) at Sa varied from 0.02 to 0.04 for sin^ 2^^23 = 0.90 
(compared to 0.03 when 6*23 = 7r/4), depending on whether 5mli is positive or negative, and 
whether ^23 < 7r/4 or ^23 > 7r/4. The sgn{5mli) sensitivities for different possibilities are 
shown in Table II. 



TABLE 11. sin^ 2^13 reach for determining the sign of Sm^i at 3a using u and i> data from JHF 
at 295 km and NuMI at L2 = 890 km with 9oa = 0.74°, for different combinations of sgn{Sm^i) 
and ^23, where sin^ 2^23 = 0.9. 





sin^ 2^13 reach for sj 


?n((5mii) 






sin 6*23 = 0.585 




sin 6*23 = 0.811 


+ 


0.04 




0.02 




0.03 




0.03 



C. Establishing the existence of CP violation 

An important goal of long-baseline experiments is to determine whether or not CP is 
violated in the leptonic sector. In order to unambiguously estabhsh the existence of CP 
violation, one must be able to differentiate between (5,^13) and all possible {S',9[^), where 
6' = 0° or 180° and 6'i^ can take on any value. For our CP violation analysis we vary 
sin^ 26*^3 in 2% increments, as was done in the previous section when testing the sgn{6m1^) 
sensitivity. 

Figure 4 shows contours of sin^ 2^13 reach for distinguishing 5 — 90° from the CP con- 
serving values S — 0° and 180° at 3a (with the same sgn{Smli)) , plotted in the {9oa,L2) 
plane, assuming u and u data at both Li and L2 are combined. The CP reach in sin^ 2^13 
can go as low as 0.01 for Oqa — 0.5 to 0.9°. Results for 5 = 270° are similar to those for 
S = 90°. 

Figure 5 shows the minimum value of sin^ 2^13 for which 5 can be distinguished from 

all CP conserving parameter sets with 5 = 0° and 180°, including those with the opposite 
sgn{6m1i) , at the 3o" level when Oqa = 0.74° and L2 = 890 km, for several different values of 
6mli. Figure 5a shows the reaches if data from JHF and NuMl are combined, while Fig. 5b 
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shows the reaches if data from NuMI only are used. For most vahies of S, when Srri^i is 
higher the CP effect is increased, and hence CP violation can be detected for smaller values 
of 6*13. However, there is a possibility that a CPV solution with one sgn((5m|]^) may not be 
as easily distinguishable from a CPC solution with the opposite sgn(5m3^); this occurs, e.g., 
in Fig. 5a for Sml^ = 1 x 10"^ eV^, where the predictions for {5 = 45° and 135°, 5mli > 0) 
are close to those for {S = 0° and 180°, 5m|^ < 0); in this case the CP reach for those values 
of 5 is about the same for Sml^ = 1 x lO^'^ cV^ and 6mli = 5 x 10~^ eV^. 

We note that if data from only JHF are used (and assuming sin^ 26*13 ^ 0.1) no value of 
the CP phase can be distinguished at 3a from the CP conserving solutions when 6ml^ < 
8 X 10~^ eV^, principally because the intrinsic CP violation due to S and the CP violation 
due to matter have similar magnitudes and it is hard to disentangle the two effects. For 
larger values of 5m^i, the intrinsic CP effects are larger and CP violation can be established; 
e.g., if 5m|i = 1 x 10"^ (2 x 10"^) eV^, maximal CP violation (5 = 90° or 270°) can be 
distinguished from CP conservation at 3(7 for sin^2^i3 > 0.006 (0.001). Therefore, when 
5m|i = 1 X 10~^ eV^, most of the CP sensitivity of the combined JHF plus NuMI data 
results from the JHF data; for Srn^i — 2x 10~^ eV^ the two experiments contribute about 
equally to the CP sensitivity. 

The boxes in Figs. 2 and 4 indicate the values of L2 and 60 a for which the cos 5 terms 
in the average probabilities vanish for the second experiment. As indicated in the figures, 
these detector positions are good for both distinguishing sgn(5m3^) (see Fig. 2) and for 
establishing the existence of CP violation (see Fig. 4), especially for larger values of L2. 
A good compromise occurs at Oqa — 0.74° with L2 — 890 km. In Ref. [15] it was shown 
that similar values for Oqa and L2 using the NuMI off-axis beam gave a favorable figure-of- 
merit for the signal to background ratio; our analysis shows that such an off- axis angle and 
baseline is also very good for distinguishing s^{5m\-^) and establishing CP violation, when 
combined with superbeam data at Li = 295 km. 

D. Resolving the (5, it — 5) ambiguity 

If L2 — 890 km and Oqa — 0.74° are chosen for the location of the second experiment, as 
suggested in the previous section, then both the first and second experiments are effectively 
measuring sin 5, and it is impossible to resolve the {S,tt — S) ambiguity. Different values of 
L2 and 9oA would be needed to distinguish 5 from tt — S. 

Figure 6 shows contours (in the space of possible L2 and Oqa) for the minimum value of 
sin^ 2^^i3 needed to distinguish 6 = 0° from 6 = 180° at the 2a level using u and P data from 
Li and L2 (it is not possible to distinguish 5 = 0° from 5 = 180° at the 3a level for any value 
of sin^ 2^13 < 0.1). Two choices are possible: one with Oqa ^ 0.3-0.5° and L2 ~ 650-775 km, 
and another near Oqa — 1-0° with L2 — 950 km. The former choice does not do well in 
distingishing sgn[Sm1i), while the latter choice is nearly optimal for sgn((5m|^) sensitivity 
but significantly worse for CP violation sensitivity. Thus the ability to also resolve the 
{d,7T — 6) ambiguity is rather poor, and comes at the expense of CPV sensitivity. 
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E. Resolving the (^23,7r/2 — ^23) ambiguity 



If 6*23 7^ 7'"/4, there is an additional ambiguity between 6*23 and 7r/2 — 6*23. This ambiguity 
gives two solutions for sin^2^^i3 whose ratio differs by a factor of approximately tan^^23, 
which can be as large as 2 if sin^ 2^23 = 0.9 [11]. Assuming Li — 295 km for the first 
experiment, we could not find any experimental configuration of L2 and Oqa for the second 
experiment that could resolve the {623.11/2 — 623) ambiguity for sin^ 2^^i3 < 0.1 at even the 
la level for the entire range of detector sizes and source powers listed in Table 1. Therefore 
we conclude that superbeams are not effective at resolving the [9 22,, 7^/2 — 623) ambiguity 
using i/g and Pg appearance data. Since the approximate oscillation probability for Ug — > 1/^ 
is given by the interchanges sin ^23 ^ cos023 and S —S in the expression for the — > i/g 
probability, a neutrino factory combined with detectors having tau neutrino detection capa- 
bility provides a means for resolving the (6*23,71/2 — 6*23) ambiguity [11]. Another possibility 
is to measure survival of Pg's from a reactor, which to leading order is sensitive to sin^ 2^13 
but not ^23 [33,34]. 

F. Dependence on jf^ml^l 

The foregoing analysis assumed l^^sil = 3 x 10~^ eV^. If the true value differs from 
this, then to sit on the peak (where the cos S terms vanish) requires tuning the beam energy 
and baseline according to the measured value of JHF has the capability of varying 

the average E^, from 0.4 GeV to 1.0 GeV, which would correspond to reahzing the peak 
condition for |5m|i| = 1.6-4.0 x 10~^ eV^ [30]. In principle, NuMI can vary both L2 and 
60A to be on the peak. If < 3 x 10^^ eV^, then the best sensitivity to sgn{6m1j^) is 

obtained for larger Oqa and longer distances (the larger angle makes smaller while the 
longer distance enhances the matter effect), and the sensitivity is reduced (since the matter 
effect is smaller for smaller Smli). The CP violation sensitivity is also reduced, although 
not as significantly. For larger values of |5m|i| the sensitivity to sgn(5m|j) is better, with 
CP violation sensitivity about the same. 

The tuning of the experiments to the peak (where the cos 5 terms in the average prob- 
abilities vanish) requires knowledge of |5m3^| before the experimental design is finalized. 
The values of |5m|i| and ^23 will be well-measured in the survival channel — > i/^ mea- 
surements that would run somewhat before or concurrently with the appearance measure- 
ments being discussed here, but of course this information may not be available when the 
configurations for the off-axis experiments are chosen. If l^mg^l is known to 10% at 3a 
(the expected sensitivity of MINOS), then the sensitivities to sgn(5m3^) and CP violation 
are not greatly affected by baselines that are slightly off-peak. If the baselines and neu- 
trino energies for the superbeam experiments must be chosen before a 10% measurement 
of can be made, a loss of sensitivity to sgn{dm1i) could result by not being on the 

peak. For example, if the experiments are designed for l^mg^l = 3 x 10~^ cV^ but in fact 
l^mg^l = 2.5 X 10~^ eV^, the 3a sgn{Smlj^) reach is less (sin^ 26*13 = 0.04, compared to 0.03 
for |(5m|i| = 3x 10~^ eV^). If \Sml^ \ is actually 2 x 10"^ eV^, the 3a sgn((5m|i) reach extends 
only down to sin^ 2^13 ~ 0.075, just a little below the CHOOZ bound. 

Since the sgn{Sm^i) determination has the worst reach in sin^2^i3 (compared to the 
discovery reach and the CPV sensitivity), and since not knowing sgn(5m3^) can induce a 
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CPV/CPC ambiguity, the measurement of sgn(5m|i) is crucial. If |5to3i| is not known 
precisely, then the exact peak position is not known, and an off-axis angle and baseline 
should be chosen that will give a reasonable reach for sgn(5m|i) over as much of the allowed 
range of |(^m3^| as possible. For example, Oqa — 0.85°-0.90° and L ~ 930 km gives a 
sgn(5m3i) reach that is fairly good for the range \5m\-^\ = 2 x 10~^ eV^ to 4 x 10~^ eV^. 
The reach for sgn(5m|]^) is farthest from optimal at the extremes (sin^ 2^13 = 0.06 versus 
the best reach of 0.05 when = 2 x 10^^ eV^ and 0.03 versus the best reach of 0.02 

when |5m|]^| = 4 x 10^^ eV^). But the CPV reach remains at least as good as the sgn(5m|]^) 
reach for this range of l^wigj^l. 

IV. SUMMARY 
We summarize the important points of our paper as follows: 

(i) Two superbeam experiments at different baselines, each measuring z/g and 

Vf. appearance, are significantly better at resolving the sgn(5m|i) ambiguity than one 
experiment alone. Using beams from a 4.0 MW JHF with a 22.5 kt detector 2° off axis 
at 295 km and a 1.6 MW NuMl with a 20 kt detector 0.7-1.0° off axis at 875-950 km, 
sgn(5m|]^) can be determined for sin^ 2^^i3 > 0.03 if ^rn^^ = 3 x 10~^ eV^. Sensitivities 
for other beam powers and detector sizes are given in Table I. 

(ii) For the most favorable cases, a higher value for the solar oscillation scale ^m^i does 
not greatly change the sensitivity to sgn(5m3-|^) when v and v data from two different 
baselines are combined (unlike the single baseline case, where the ability to determine 
sgn(5m|i) is significantly worse for bm^^x ^ 5 x 10~^ eV^). 

(iii) Running both experiments at the oscillation peaks, such that the cos b terms in the 
average probabilities vanish, provides good sensitivity to both ^^{bm\^ and to CP 
violation. On the other hand, the ability to resolve the (5, tt — 5) ambiguity is lost, 
and the (^23,7r/2 — ^23) ambiguity is not resolved for any experimental arrangement 
considered. However, the (5, tt — 5) and (^23, '?'"/2 — ^23) ambiguities do not substantially 
affect the ability to determine whether or not CP is violated (although the latter 
ambiguity could affect the inferred value of 6'i3 by as much as a factor of 2). 

(iv) Since running at or near the oscillation peaks is favorable, knowledge of |5w^3i| to 
about 10% (from MINOS) before these experiments are run would be advantageous. 
If |5m3i| is not known that precisely in advance, then the detector off- axis angle and 
basehne can still be chosen to give fairly good (though not optimal) sensitivities to 
sgn(5m|i) and CP violation. 

We conclude that superbeam experiments at different baselines may greatly improve the 
prospects for determining the neutrino mass ordering in the three-neutrino model. Since 
a good compromise between determining %'ga.{bm^^ and establishing the existence of CP 
violation is obtained when both experiments are tuned so that the cos 8 terms in the average 
probabilities approximately vanish, knowledge of would be helpful for the optimal 

design for the experiments. 
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FIG. 1. Contours of (a) best-case (when Sm^^ > 0), and (b) worst-case (when Sm'^i < 0), 
sin^ 2^13 3a discovery reach in the {O0A7L2) plane, for the v channel at NuMI, where Oqa is the 
off-axis angle and L2 is the baseline of the NuMI detector. For the other neutrino parameters we 
assume \6mli\ = 3 x 10"^ eV^, 6*23 = vr/4, 5mli = 5 x 10"^ eV^, and sin^ 26*12 = 0.8. The boxes 
indicate detector positions for which the cos 6 terms in the average oscillation probabilities vanish. 
For the i> channel the results are similar, except that the best case occurs for Sm^i < and the 
worst case for Sm^i > 0. 



14 



1000 



400 



sin^(2 9,3) = 0.030 

0.040 

0.060 
0.080 








































□ 












□ 

□ 


□ 


□ 



0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

eoA(°) 

FIG. 2. Contours of sin^ 2^13 reach for resolving the sign of drn^i at the 3a level in the {OqAj L2) 
plane when data from JHF and NuMI are used. The JHF detector is assumed to have baseline 
Li = 295 km. Other parameters and notation are the same as in Fig. 1. 
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FIG. 3. Minimum value of sin^ 2^13 for which sgn{6m'^i) may be determined at 3cr, assuming 
the true solution has Sm'^i > 0, using v and P data from (a) JHF with Li = 295 km and NuMI 
with L2 = 890 km, and (b) only NuMI with L2 = 890 km, for several values of Sm2i (in eV^). The 
off-axis angle for the NuMI detector is 60 a = 0.74°. Other parameters are the same as in Fig. 1. 
Results for Sni^i < are approximately given by reflecting the curves about S = 180°. 
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FIG. 4. Contours of sin^ 2^13 reach for distinguishing S = 90° from the CP conserving values 
5 = 0° and 180° at 3a (for the same sgn{Sm^i)), plotted in the (^oA)-^2) plane, when data from 
JHF and NuMI are combined. Other parameters and notation are the same as in Fig. 1. Results 
for 6 = 270° are similar to those for 6 = 90°. 
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FIG. 5. Minimum value of sin^ 20i3 versus CP phase 5 for which 5 can be distinguished from 
the CP conserving values (5 = 0° and 180° (with either sign of Smi^i) at the 3cj level when (a) data 
from JHF and NuMI are combined, and (b) data from NuMI only are used. The baseline for JHF 
is Li = 295 km, while for NuMI L2 = 890 km and Oqa = 0.74°. The curves are plotted for several 
values of the solar mass scale ^m^i (in eV^). Other parameters are the same as in Fig. 1. 
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FIG. 6. Contours of sin^ 2^13 reach for distinguishing (5 = 0° from 8 = 180° (the ((5, vr — 5) 
ambiguity) at the 2a level, plotted in the {9oa,L2) plane, when data from JHF and NuMI are 
combined. Other parameters and notation are the same as in Fig. 1. 
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